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Difference of degeneracy of the low-spin (LS) and high-spin (HS) states causes interesting 
entropy effects on spin-crossover phase transitions and charge transfer phase transitions in 
materials composed of the spin-crossover atoms. Mechanisms of the spin-crossover (SC) 
phase transitions have been studied by using Wajnflasz model, where the degeneracy of the 
spin states (HS or LS) is taken into account and cooperative natures of the spin-crossover 
phase transitions have been well described. Recently, a charge transfer (CT) phase transition 
due to electron hopping between LS and HS sites has been studied by using a generalized 
Wajnflasz model. In the both systems of SC and CT, the systems have a high temperature 
structure (HT) and a low temperature structure (LT), and the change between them can be a 
smooth crossover or a discontinuous first order phase transition depending on the parameters 
of the systems. Although apparently the standard SC system and the CT system are very 
different, it is shown that both models are equivalent under a certain transformation of 
variables. In both systems, the structure of metastable state at low temperatures is a matter 
of interest. We study temperature dependence of fraction of HT systematically in a unified 
model, and find several structures of equilibrium and metastable states of the model as 
functions of system parameters. In particular, we find a reentrant type metastable branch 
of HT in a low temperature region, which would play an important role to study the photo- 
irradiated processes of related materials. 



§1. Introduction 



In the so-called spin-crossover atoms, e.g., Fe and Co, spin states of the atoms 
can be changed between the high-spin (HS) state and the low-spin (LS) state by 
small perturbations, such as changes of temperature and pressure and also photo- 
irradiation, etc. This change of spin state is called spin-crossover (SC) transi- 
tion. 1 ) -41 ) The HS state has a large number of degeneracy and it is favorable at 
high temperatures, while the LS state has low energy which is favorable at low tem- 
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peratures. The competition between the effects of degeneracy (the entropy) and of 
the energy provides interesting changes of structure of material. 

Prom a view point of adiabatic energy levels of local structure of atoms, the 
population dynamics among states (HS, LS, and other excited states) under photo- 
irradiation has been studied by a kind of rate equation, which explains mechanism 
of transition among the LS and HS states. 15 )' 16 ) It has been also pointed out that 
cooperative interactions are important for SC transitions. With an aid of interaction 
between atoms, this transition can be either smooth crossover or discontinuous first 
order phase transition depending on the system parameters. 1 )~ 4 )' 21 )- 34 )> 4 °)> 41 ) This 
situation is well described by the Wajnflasz model. 40 )' 41 ) Control between the HS 
and LS states has been realized by photo-irradiation as the light-induced excited 
spin state trapping (LIES ST), 1 3)-20),35)-39) anc j structure of the metastable ordered 
state of the systems has become an important topic. 20 )' 42 )~ 45 ) 

Prussian blue analogues (PBAs), classified as molecule-based magnets, 46 ) -53 ) 
show various photomagnetic phenomena. 54 )~ 62 ) Two-way photo-switching between 
magnetic and nonmagnetic states was observed in a Co-Fe Prussian blue analogue, 
Ko.2Coi.4Fe(CN)66.9H20, 55 )' 56 ) which has attract much interest. This photomag- 
netism is based on the charge-transfer-induced spin transition (CTIST) between two 
phases, 58 ) e.g., one is the high temperature (HT) phase consisting mainly of the 
Co n -Fe m where Co 11 is in the high spin state (HS, 5 = 3/2) and Fe m is in the low 
spin state (LS, 5 = 1/2), and the other is the low temperature (LT) phase consisting 
mainly of the Co ni -Fe n where Co 111 is in the low spin state (LS, 5 = 0) and Fe 11 is 
also in the low spin state (LS, 5 = 0). In this material, Co and Fe ions are antiferro- 
magnetically coupled. Therefore, precisely speaking, it is a ferrimagnetic material. 
However, the spin-crossover transition occurs at the Co site, and if we take a unit 
of the pair atoms (Fe and Co), we may regard the transition as a generalized SC 
transition. 

On the other hand, another type of phase transition due to the charge trans- 
fer (CT) was discovered in a mixed-valence iron complex (n-C3H7)4N[Fe n Fe m (dto)3] 
(dto = C202S2). 63 )~ 65 ) In this material, the charge transfer causes a change of degen- 
eracy of atomic states, and a phase transition occurs between the high temperature 
(HT) and low temperature (LT) structures. In this system, the spin transition on 
each atom does not occur. However, the entropy effect has an important role for 
the structure change of electron configuration and nature of phase transition was 
explained by using a generalized Wajnflasz model. 65 ) There, it was found that a 
metastable branch of HT exists at all the temperatures below the critical tempera- 
ture. This behavior is qualitatively different from that in the standard model of SC 
transition. That is, because experimentally the fraction of HS shows a hysteresis loop 
in SC transitions, for SC transitions we have adopted a model where the metastable 
branch of HS terminates at a spinodal point, and below which the HS is unstable. 
In this model, it is difficult to explain the metastable behavior of photo-induced 
HS state at low temperatures. In order to overcome this difficulty, a dynamical 
effect has been introduced by adopting the so-called Arrhenius dynamics, 35 ) which 
has explained successfully the fact that the HS state remains for a long time at low 
temperatures although it is unstable in a mean-field free energy. 37 ) -39 ) On the other 
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hand, in the model of CT, the HT state is metastable at all the temperatures. Now, 
it is an interesting problem to study the relation between the both models. 

In this paper, we will show that the both models are equivalent under a certain 
transformation of variable, and the unified model shows several qualitatively different 
structures of equilibrium and metastable states as a function of system parameters. 
In particular, a reentrant type metastable branch of HT in a low temperature region 
is discovered. The existence of this low temperature metastable structure would play 
an important role to study the photo-irradiated processes of related materials. Here, 
we mainly study the transition of spin states, e.g. HS 44> LS. However, if we consider 
magnetic interaction between the spins, the structure of phase transition is modified. 
The magnetic effect has been studied for Co-Fe PBA. 38 )' 39 )' 43 )' 44 ) The magnetic 
phase transition in the CT system has been also studied. 65 )~ 67 ) The existence of 
the metastable states causes new variety of magnetic phase transition, which will be 
briefly discussed. 

§2. Unified model 

2.1. Typical spin-crossover phase transition with a single hysteresis 

In order to describe the spin-crossover transition, the so-called Wajnflasz model 
has been adopted, 40 -*' 41 -* where the HS state is represented by s = 1, and the LS 
state is represented by s = — 1. It should be noted that the states Sj = ±1 have 
different degeneracy. Let the number of states of the HS state be bh, and that of the 
LS state til- Let us review the Wajnflasz model briefly. 35 ) The interaction among 
atoms is originated in the elastic interaction in the atomic structure. 25 ) However, 
here it is simply modeled by a nearest-neighbor interaction between sites, and the 
Hamitonian has the following form 

<ij> i 

where J denotes the interaction between states of neighboring atoms, and D denotes 
the energy difference between HS and LS. Here we consider the case where D > 0, 
i.e., the LS state is energetically favorable. The partition function is given by 

z = Xre -/3H w= J- £ 'expipj^SiSj-pDY^sA, (2-2) 

si=±l sjv=±1 V <ij> i I 

where Yl' denotes that the summation is carried out over the degenerate states, i.e., 
summation over nn states of Sj = 1 and states of Sj = —1. Here, f3 = 1/ksT, and 
N is the number of sites. This form is expressed by a non-degenerate Ising variable 
<jj = ±1 as 



Z= £ • • • £ exp i(3J J2 Wj ~ PD Y °i II ("h^.i + niAr.-i) • (2-3) 

<ri=±l <jjv=±l V <ij> i ) i 
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We may rewrite the second factor as 



lo-ln 2iL 

^H^i + = ^nnn L e 2 ™l . (2-4) 



Now we have 



Z = (n H n L ) 2 ^ • • • ^ exp Vi<?j - (d - In 5 j ^ a. 

CT1=±1 crjv=±l \ \ <ij> ^ ~ i 



(2-5) 



and 



9 = — • (2-6) 
til 

Therefore the model (|2-1[) can be expressed by an effective Hamiltonian with a tem- 
perature dependent field 

n = - J E + ( D - ln9 ) ^ a " (2 ' 7) 

This form of Hamiltonian is called Wajnflasz model. 

Here, it should be noted on a characteristic property of this model. At high tem- 
peratures, k^T > 2D/\ag, the term of the effective field (-ff e ff = — D + k^T In g/2) is 
positive, and thus the spins has a positive expectation value, (<7j) > 0. On the other 
hand, at low temperatures, k-^T < 2D/lng, (<jj) < 0. Let us define the marginal 
temperature To to be 

k B To = (2-8) 
In g 

where (<7j) = 0. If To is larger than the critical temperature T\c of the corresponding 
Ising model 

<ij> 

that is, 

T > Tic, (2-10) 
then the change of the magnetization is smooth. On the other hand, if 

T < T IC , (2-11) 

then a discontinuous change occurs, which means the first order phase transition 
takes place. 

A typical example of the hysteresis obtained by a kind of mean field analysis is 
schematically depicted in Fig. ^ So far, the experimental results of the first order 
phase transition have been explained by this type of temperature dependence of order 
parameter /hs (= ((fi) + which is the fraction of the atoms in the HS state. In 
this figure, the HS state becomes unstable at a temperature T 1 denoted by a dotted 
line. Below this temperature, the HS state is unstable. Thus, it seems difficult 
to have a long-lived HS state at low temperatures. In Co-Fe PBA, however, the 
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Fig. 1. Schematic temperature dependence of a solution of the self-consistent equation of a mean- 
field theory for the high spin fraction q of models of SC transition with a hysteresis loop. The 
bold dotted curves denote the hysteresis region. The thin dotted line denotes the first order 
phase transition temperature where q = 0.5. 



magnetic state exists for a long time after pumping by photo-irradiation. In order to 
explain this long-lived state we have introduced a mechanism of very slow dynamics, 
i.e., the so-called Arrhenius dynamics, 35 )' 39 ) where we assume an Arrhenius type 
relaxation time: r oc exp(— E^/k^T) with an appropriate value of Eq. 

2.2. Metastable structure of a charge transfer phase transition 

Contrary to the above mentioned behavior, we have found a case where the 
static metastability exists in the study of the charge transfer transition in material 
(nC3H 7 )4N[Fe n Fe m (dto)3] (dto=C 2 2 S2). 63) ~ 65) Here, let us consider this charge 
transfer transition. This material consists of bipartite lattice. At one of the sublat- 
tices, Fe is surrounded by sulphur atoms where Fe is always in the low spin state 
(LS), and at the other sublattice Fe is surrounded by oxygen atoms where Fe is al- 
ways in the high spin state (HS). We call the former 'A-site' and the latter 'B-site'. 
Let the both sublattices have N sites. Fe 11 has one more electron than Fe 111 , and the 
difference between them can be expressed by the number of additional electron rtj 
which is 1 in Fe 11 and in Fe 111 . Here, we consider the degeneracy of spin degree of 
freedom. 65 ^ At A-site, Fe 11 is in the low spin state and S = as depicted in Fig. [2] 
and thus the degeneracy is 1, and the spin of Fe 111 is S = 1/2 and thus its degeneracy 
is 2. Similarly at B-site, 5 = 2 for Fe 11 and S = 5/2 for Fe 111 . In terms of {ni\, the 
degeneracy is given as the following: 



site 


/; 




spin 5* 


degeneracy 


A 


n — 


1 





UAl — 1 
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S=2 S=5/2 S=0 S=1/2 





^ ^ ^ 

^ ^ 

Fe(ll) Fe(lll) Fe(ll) Fe(lll) 

(a) High spin (HS) states (b) Low spin (LS) states 

Fig. 2. Electron states of Fe 11 and Fe m (a) in the high spin (HS) state and (b) in the low spin 
(LS)state 

Because the numbers of sites of Fe 11 and Fe m are the same, the total number of the 
additional electrons is N 

~^m = N. (2-12) 



In this material, electrons transfer between A and B sites. We have introduced the 
following Hamiltonian for this system 

H ct = £ Yl nifij + A^m, (2-13) 

<ij> i€B 

where e represents the interaction between the electrons at nearest neighbor sites. 
A{> 0) denotes a kind of local on-site energy. Because of this term, energetically 
electrons tend to stay at A-sites. If all the electrons are on A-sites, then the system 
does not cost energy due to A. Therefore, we define this state as the perfect LT 
state. Because of the difference of the degeneracy listed above, at high temperature, 
electrons tend to stay at B-sites. We express the degree of high temperature (HT) 
state by the number of electrons at B-sites 

<n B > = (2-14) 

In the above choice of the degeneracy, (ne) is more than half at high temperatures. 
In the present model, no spin-crossover transition occurs by the electron transfer. 
However, the degeneracy of spins still plays an important role, and this model ex- 
hibits a first order phase transition. 65 ) In the mean field theory, the self-consistent 
equation for q = (<7i 6 b) is given by 



tanh 



- j + fcBrin(5/3))j 



(2-15) 



where z is the number of the nearest neighbor sites and is 3 for the Honeycomb 
lattice. We depict the typical temperature dependence of solution of self-consistent 
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equation for (rae) = (<Z+ l)/2 in Figs.EJfor A = 1 and 10. There, we fix e to be 3. In 
Fig-Etb), there are three solutions for q at low temperatures. There, the largest and 
smallest ones give stable solutions and the second one gives an unstable solution. 
By comparing the free energies corresponding to them, we find the first order phase 
transition point, which is To. In figures, the temperature where (h-b) = 0.5 is shown 
by a dotted lines. This temperature has the same physical meaning as to To defined 
in (|2-8|) as will be shown in the following subsections. 

Here we find again that there are cases where the change is smooth and discon- 
tinuous. However, in this charge transfer transition, the metastable branch of HT 
state remains until T = 0, i.e., HT state is metastable even at very low temperatures. 




(a) (b) 

Fig. 3. Temperature dependence of fraction of the HT state (tib) for systems with e = 3. (a) 
A = 10 and (b) A = 1. To is the temperature where (jib) = 0.5 and also it is the critical 
temperature in the case (b). 



2.3. Comparison 

In the two observations in the above SC and CT cases, there are qualitatively 
different temperature dependences of the free energy structure. Thus, it would be 
interesting to compare these two cases, and we will attempt to express this model 
(|2~T3^) in the form of Wajnflasz model 

First, we transform the variable n,(= or 1) to Sj(= —1 or 1): 

ni = *±i. (2-16) 
Here the constraint = N) becomes 



i 



(2-17) 
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The Hamiltonian becomes 



n = i E (« + + x ) + 1 E( s * + !)• ( 2 - 18 ) 

<ij> ieB 



Making use of the constraint (|2-17j) we can rewrite it as 

n = 1 E & + ^ + x ) + j ( E(^ + !) - + !) J • ( 2 - 19 ) 

<ij> ViGB iGA / 

In order to have a ferromagnetic model, we perform a local gauge transformation: 
s — > — e, Sj — ► — 5j at A-site and s, — > Si at B-site. Then we have 



<2jr" 



E^ + j(E^ + E^)- ( 2 - 2 °) 

<M> VtGA iGB / 



Here we again use non-degenerate Ising variable o~j. It should be noted that the 
numbers of degeneracy at A- and B-sites are different. However, using relations 
similar to l)2-4|) . the numbers of degeneracy (nAi, hao, ^bi and neo) can be taken 
into account as 

n = 'I + (t + ^ BTln ^) X> + (t - ^ BTln S) i> 

" (2-21) 

Using the constraint 

E^~E^ = > ( 2-22 ) 

iGA iGB 



E^+^ + I^Tln^) J>, (2-23) 



we finally have 

<«i> 

where oa = ^ai/^ao an d 5b = ^bi/^bo- Now, we have the same form as that of 
Wajnflasz model (|2-7|) . Therefore, it is proved that the model (|2- 13|) is equivalent to 
the model (j2~T|) . 

2.4. To /or i/ie CT transition 

In this formalism, we can make use of the relation Q2-8|) to distinguish the smooth 
and discontinuous changes. Here, the transition temperature Tq is given by 

A 
In" 



fc B7o = (2-24) 



9A 

Substituting the numbers of degeneracy, we have 



* bT » = d« = tali * L958A (2 - 25> 

111 172 ' 
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In this subsection, we set ks = 1. 

In the case A = 10, the condition (12-10(1 is satisfied. That is, in the honeycomb 
lattice, the critical temperature of the ferromagnetic Ising model is given by fceTci = 
3(e/4) in the mean-field theory, and therefore, for e = 3, 

T C i = 0.75 x 3 < T = 1.958 x 10, (2-26) 

and thus (tib) changes smoothly, which is depicted in Fig. Ota). On the other hand, 
for A = l, 

T C i = 0.75 x 3 > T = 1.958 x 1, (2-27) 
and the transition is of the first order, which is depicted in Fig. Efb). 

2.5. Monte Carlo study 

It should be noted that, although in the above we study the model in the mean- 
field approximation, the criterion whether the change of (tib) at To is smooth or 
discontinuous, i.e., the relation between To and Tci, also holds in the exact treatment. 
In the above treatment, the value of To is given in the mean-field approximation 
because we studied the model in the mean-field approximation. However, in the 
exact treatment, To is given by the critical temperature of the ferromagnetic Ising 
model: 

T C i = 2 —^-- t ^ 1-52 x (i) . (2-28) 

ln(2 + ^3)4 W v ; 

In Fig^J we depict the temperature dependence of (tib) studied by a Monte Carlo 
(MC) method for e = 3 and 6 with A = 1. Because we study the exact short range 
model in MC, we find exact properties in MC although it is numerically done. For 

s = 3, 

T C i = 1.52 x (3/4) < T = 1.958, (2-29) 

and then we indeed find that (tib) changes smoothly as plotted by the blue dots. 
For discontinuous change of (tib), the relation To < Tci, i.e., 

1.958 < 1.52 x (e/4), (2-30) 

must be held. Therefore, in MC, e must be larger than 5.15(= 1.958 x 4/1.52) for the 
discontinuous change. The temperature dependence of (tib) for e = 6 is also plotted 
in Fig. 0]with black dots, where we find a discontinuous jump. Here, we changed the 
value of e instead of A to have the same value of To for the convenience to plot the 
data in a figure. 

Monte Carlo simulations were performed on the honeycomb lattice of the size 
32 x 32. The simulation started at T = 0.125 where (tib) — in the equilibrium state. 
The temperature was increased up to T = 4. At each temperature, we performed 
20000 Monte Carlo steps (MCS) for a transient process and then took data in the 
subsequence 100000MCS. Fluctuation of the data is smaller than the size of the dots. 

We find that (tib) jumps at T = 2.12 although the exact critical point is To = 
1.958 which is indicated by the dotted line where (tib) = 0.5. This shift of the jump 
is due to the hysteresis phenomenon associated with the first order phase transition. 
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After T = 4, the temperature was reduced. We find almost the same values 
of (n B ) as those in the heating process above T = 2.12. However, it keeps a large 
value below this temperature, and stays in metastable HT state even below To. This 
metastability corresponds to the metastable solution found in the mean field theory. 
In the case e = 3, no hysteresis is found. 




T 



Fig. 4. Temperature dependence of fraction of the HT state (ne) for e = 3 (blue dots) and 6 (black 
dots) studied by Monte Carlo simulations. Here, A = 1. 

In this way, we can know various properties of (n B ) from the knowledge on the 
corresponding Ising model. For example, at the critical value of e, (n B ) shows a 
second order phase transition where 

(n B )-^oc|r-r |7, 6 = 15 (2-31) 

and above the critical value of e, («b) changes discontinuously. The jump of (n B ) is 
given by 

(n B (T + 0)-n B (T -0)> oc \e - e c f , /3 = ~, (2-32) 

o 

in the exact treatment and also in the Monte Carlo simulation. In the mean-field 
theory, on the other hand, the critical exponents are 5 = 3 and (3 = 1/2. In Fig. 
the temperature dependences of (n B ) at the critical value of e are plotted for both 
cases studied in the mean-field theory (by bold curve) and in MC (by circles). Here 
it should be noted that (n B ) shows a non-monotonic temperature dependence even 
in the case of the second order phase transition in MC. 

§3. Classification of temperature dependence of the high-spin fraction 

In the previous section, we have found the equivalence of the model of the SC 
transition and the model of the CT transition. In the SC case, the HS state becomes 
unstable at a temperature T' below the hysteresis loop. On the other hand, in the 



Structure of Metastable States in Phase Transitions 



11 




Fig. 5. Temperature dependence of fraction of the HT state (tib) at the critical value of e studied 
by the mean-field theory ec = 3/(4 In 10/3) = 2.61 • • • and by a Monte Carlo simulation ec = 
1/(4 In 10/3) = 5.18- • • . Here, A = 1. 



CT case the HT state remains metastable below To until T = 0. These two observa- 
tions suggest that there are qualitatively different temperature dependences of the 
free energy structure. By studying these structures, we can reach a comprehensive 
understanding of phase transitions of systems of spin-crossover atoms. 

Here we study the model of the type of 1)2 -7|) with the variable cr, = ±1 

k B T 



H 



<*3> 



D 



<7,;. 



(3-1) 



As has been mentioned in the previous section, the temperature T$ where (o"j) = 
is given by (|2-8j) . Here, we define Dq as that the system exhibits a first-order phase 
transition when D < Dq: 

Dc = ^m g . (3-2) 

Next, we consider the condition for the metastability of the HT state at the ground 
state. By flipping one spin in the configuration where all the spins are +1, the system 
gains the crystal field —2D, while it loses the exchange energy 2zJ where z is the 
number of nearest neighbors. Thus, if D is smaller than the critical value: 



Dqg = zJ, 

the all up state is metastable at T = 0. 

Let us study temperature dependence of solution(s) 

q(t) = (en) 

of mean-field self-consistent equation: 



tanh 



P Jzq 



D + ^lng 



(3-3) 



(3-4) 



(3-5) 
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T/D 



Fig. 6. Temperature dependences of q(T) are depicted for D — 10.0 (thin solid curve), D — 7.2 
(dashed curve), and D — 6.3 (bold solid and dashed curve). Here, J = 1 and a = 1.2 (g = e 2 ' 4 ). 



F(q,T) 



F(q,T) 





(a) 



(b) 



Fig. 7. Temperature dependence of the free energy of functions of q are depicted for (a) D = 10.0, 
and (b) D — 6.3. The temperature changes by 1 from T = 15 (red) to 1 (blue). Here J = 1 and 
a = 1.2 (g — e 2,4 ). The open circles indicate extreme points of the free energy. 



The free energy for the solution is given by 



F(q) = -zJq 2 -k B T In 



, . Jzq — D 

2 cosh | h a 

k B T 



(3-6) 
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Here, we adopt the following parameters 



J = 1 



z 



= 6, and a= -lag = 1.2 



(3-7) 



and here we take J as a unit of the energy. In the followings, we study the temper- 
ature dependence of q(T) and its dependence on D. 

3.1. Large D behavior 

In Fig. |SJ we show the temperature dependence of q for D = 10 by a thin 
solid curve where we find a smooth change of q, and corresponding temperature 
dependence of the free energy F(q) — F(0) is depicted in Fig. Ufa). In the present case 
Dq is 7.2, below which the system shows a first-order phase transition. Temperature 
dependence of q(T) for D = 7.2 is depicted by a dotted curve. The temperature 
dependence of q for D = 6.3 is depicted by bold solid and dashed curve, which shows 
a first order phase transition. The overhanging part is plotted by a bold dotted 
line. This temperature dependence corresponds to the standard phenomena of HS- 
LS transition shown in Fig. ^ Corresponding temperature dependence of the free 
energy is depicted in Fig. E[b). 



Fig. 8. Temperature dependences of q are depicted for D = 5.2 (bold solid curve) and D = 4.0 
(thin solid curve). The dotted curves are q(T) for D = 10, 7.2, and 6.3 which are depicted in 
Fig. E| Here, J = 1 and a = 1.2 (g = e 2A ). 

In typical spin-crossover transitions, the ratio of the degeneracy g is considered 
to be rather large, and it causes the value of q in the high temperature limit 



q(T) 



-1± 



a 




q(oo) = tanh(a) 



9-1 



(3-8) 



9 + 1 



to be large. 
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F(q,T) 



F(q,T) 




1 q 



(a) 



(b) 



Fig. 9. Temperature dependence of the free energy of functions of q for (a) D — 5.2 and (b) D — 4.0. 
The temperature changes by 1 from T = 15 (red) to 1 (blue). The data for T = 0.5 is also 
plotted by a bold dotted blue curves. Here, J = 1 and a — 1.2 (g = e 2 ' 4 ). 



3.2. Hysteresis branch in low temperature region for the intermediate D 

If we take the value of D below Dqq, we find a metastable branch appears at 
low temperatures. In Fig. |SJ we show the temperature dependence of q for D = 5.2 
by a bold curve, and corresponding temperature dependence of the free energy in 
Fig. EJa) . Here we find the free energy has a local minimum of HS at low temper- 
atures. This low temperature metastable branch of HS may play an important role 
for the long-lived metastable state in Co-Fe PBA. In this case, the metastability of 
HS is intrinsic and if the initial state is set to be near the metastable point, the 
value of q is expected to move first to the metastable value and then relaxes to the 
LT state through a nucleation process. On the other hand, if there is no metastable 
branch at low temperatures as the case D = 6.3, the long-lived HS state must be 
purely due to the slow dynamics. In this case, q always moves to the 0. Checking 
this initial move of q, one may find which case is realized. It would be an interesting 
problem to determine that to which case individual materials belong. 

3.3. Small D behavior 

If we reduce D furthermore, e.g., D = 4, HT state becomes metastable at all 
temperatures as we depict in Fig. [HI by a thin curve. The free energy change of this 
case is plotted in Fig. E{b). This case corresponds to the dependence found in CT 
transition as shown in Fig. |31 For D = 4, the high temperature branch connects to 
the metastable HT branch at low temperatures, which is topologically different from 
those of large values of D. It should be noted that in this case the potential barrier 
between the metastable HT state and the stable LT state shows a non-monotonic 
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dependence, and the metastable state easily relaxes to the stable LT state via a 
nucleation process in the intermediate temperature region. The dependence of the 
barrier is plotted in Fig. I1U1 

&$E, &$E/T, q(T) 




T/D 



Fig. 10. Temperature dependence of the barrier height between the metastable HT state and the 
stable LT state for D = 4.8. The temperature dependence of q(T) is also plotted by thin dotted 
curves. The energy barrier AE is depicted by the bold dotted curve, and AE/k^T is depicted 
by the bold curve. Here, J = 1 and a — 1.2 (g = e 2 ' 4 ). 



3.4. Topological change of the structure of solutions 

So far we found two regions of D. In the large D region the solution of HT 
state at high temperatures connects to the LT state at low temperatures. On the 
other hand, in the small D region, the solution of HT state at high temperatures 
remains until T = 0. There, a LT solution exists separately. Thus, there is a critical 
value of D between these two regions. We depict q(T) at the critical value of D 
in Fig. 1111 In the mean-field theory, q(T) does not depend on the temperature at 
this marginal value of D. The values of q(T) and the marginal D are found in the 
following analysis. If the following relations 



D 

7z 



and q = tanha, 



hold, then from the self-consistent equation (|3-5j) . 

, , Jzq — D 

q = tanh( — — h a) 

k B T 



(3-9) 



(3-10) 



is always satisfied. From this relation, we find that the critical value of D is given 
by the equation: 

Dcx = z J tanha (3T1) 

and 

q = tanh a. (3-12) 
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Indeed the marginal value of D is given by 

D C x = 6 x tanh(1.2)J = 5.00193 ■■■J. 



(3-13) 



q(T) 




T/D 



Fig. 11. Temperature dependence of q(T) for D = Z?cx (solid curve). The temperature dependence 
for the values studied previously are shown by dotted lines (D = 10, 7.2, 6.3, 5.2 and 4.0). Here, 
J = 1 and a = 1.2 (g = e 2 4 ). 

3.5. Sequence of the types of structure 

Here we found three critical values of D. Finally let us point out the fact that 
the isolated low temperature metastable branch of HS exists for all the values of a. 
That is, in order to have the metastable state at low temperature, D must be smaller 
than -Dec i-e., 

D < zJ, (3-14) 

and in order to avoid that the solution of HT exists at all the temperatures, D must 
satisfy the condition 

D < D C x = -zJtanha. (3-15) 

Both conditions are compatible although the region of D is rather narrow for large a. 
In experiments, one may control the value of degeneracy by changing the pressure, 
etc. It would be very interesting to find the qualitative change of temperature 
dependence of metastable states which is found here. 

§4. Summary and Discussion 

In this paper, we studied models for phase transitions of systems consisting of 
spin-crossover atoms. We pointed out that the apparently different models for the 
spin-crossover transition and for the charge-transfer transition are equivalent, and 
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we comprehensively studied structures of stable and metastable states of the unified 
model. We found several qualitatively different ordering processes, and found critical 
values of D between different types of ordering processes. In particular, we found 
a metastable HS state at low temperatures which exists separately from the high 
temperate HS state. We find that a metastable HS state universally exists in the 
present type of models. This metastable HS state at low temperature would play an 
important in pumping process by photo-irradiation. The present study will be useful 
to classify the metastability in various photomagnetic materials where metastability 
will be investigated by photo-irradiation. 

In Co-Fe PBA and CT materials, the system shows a magnetic ordering in 
addition to the spin structure change. There, we have to consider the magnetic 
interaction. Effect of the magnetic interaction on the combined ordering process of 
fraction (q) of HS and magnetic order (m) has become an interesting problem. Here 
it should be noted that the magnetic state is always metastable but not equilibrium 
unless the magnitude of the interaction is of the order of D in the models of SC. If 
D is so large that the magnetic state appears in equilibrium, the change of HS-LS 
is smooth and the first order phase transition does not takes place. The detailed 
properties of magnetic transition will be reported elsewhere. 38 )' 39 )' 68 ) This feature 
of metastable magnetic state is compatible with the observation in Co-Fe PBA. 

In the CT system, however, magnetic ordered state appears as an equilibrium 
state at a much lower temperature than the CT transition temperature. In order to 
explain this magnetic transition, we need to consider the special property of the CT 
system. Indeed, in the CT system, magnetic moments exist on B-sites even at low 
temperatures. Thus, if some mechanism exists to connect the magnetic moments, 
the magnetic order can be formed even in the low temperature configuration. One 
of the authors has proposed a mechanism that the thermal fluctuation of A-site 
could mediate the magnetic ordering between the B-sites. 65 ) However, this mech- 
anism has been found to be difficult at least in two dimensional systems, 66 ) and 
several alternate origins of magnetic ordering are under investigation. In particular, 
it is possible that a quantum fluctuation helps the magnetic ordering which will be 
reported elsewhere. 67 ) 
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